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§ 1 . Statement of the Problem and Historical References. 

If we consider a pencil of rays issuing from a point on the axis of a symmetrical 
optical system (i.e., a system of refracting spherical surfaces, the centres of which lie 
on a straight line called the axis of the system), it is well known that, if the pencil 
be a thin one, of which the mean ray is along the axis, the first approximation to the 
emergent pencil is another punctual pencil, of which the rays pass through an image 
point, also situated on the axis. The general method of treatment of such image 
points, which are usually referred to as " geometrical" images, is due to Gauss, and is 
developed in any text book of Geometrical Optics. 

When, however, the pencil considered is one of finite aperture, the outlying rays 
do not, after emergence, pass through the Gaussian image point, nor do they have the 
inclination assigned to them by the Gaussian calculation. The emergent rays lying 
in any one axial plane touch an envelope or caustic, which has one cusp at the Gaussian 
image, with the axis as proper tangent. The intercepts of any given emergent ray 
upon the axis and the image plane, measured from the Gaussian image, are known as 
the longitudinal and transverse spherical aberrations of that ray. 

It is clear that if both these spherical aberrations, or either of them together with 
the inclination of the ray on emergence, be known for every possible position of object 
and image, and for every possible inclination of the incident ray, the whole complex 
of emergent rays lying in axial planes can be mapped out. The calculation of these 
aberrations is therefore of fundamental importance in practical optical design, where 
we do not deal with infinitely thin pencils. 

The method employed hitherto for dealing with aberrations from the mathematical 
standpoint has been to develop the sines occurring in the refraction equations at each 
spherical surface in ascending power of some argument, which may be either the 
circular measure, or the sine, or the tangent, of one of the angles concerned, and to 
calculate, by the usual methods of successive approximation, the required aberrations 
as a series of ascending powers of such argument. 
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When this is done it is found that the terms due to the first power of the argument 
lead to the Gaussian image point, so that the series begin with a term involving the 
second power of the argument in the case of the longitudinal aberration, and the 
third power in the case of the transverse aberration. These terms are the first order. 
The following terms next in sequence, which are of fourth and fifth power respectively, 
are spoken of as aberrations of the second order, and so on. 

A considerable amount of theoretical work has been done on aberrations of the 
first order by Seidel, Abbe and others, and the treatment of these is fairly well 
known. Unfortunately, it is found in practice that the first order aberrations do not 
give a sufficient approximation for the optician's requirements. In fact for a certain 
range of object and image positions, they are so badly out that they cannot be 
said to constitute an approximation at all. This fact has long been recognised by 
optical designers, whose practice is invariably to calculate, using the exact trigono- 
metrical equations which involve no approximation at all, the correct paths of a 
number of selected rays, from which they draw conclusions as to the efficiency, or 
otherwise, of the proposed system from the practical point of view. 

The trigonometrical method, however, from the designer's point of view, has the 
radical defect that, while it gives partial information about the performance of a given 
system, it gives no direct intimation of the direction in which the elimination of various 
defects is to be looked for, and it entails a long and laborious process of seeking for 
the optimum by trial and error. 

The object of the authors of the present paper has been to develop a method of 
expressing the aberrations, which, while carrying the algebraic development to a 
stage including the second order, should be free from certain grave troubles involved 
by failure of convergency, troubles which appear to have been hitherto neglected. 
In fact this method gives numerical results that, for a single lens, are considerably 
more accurate than the ordinary second order formulse. Further, these methods 
enable one to deal, in a comparatively easier form, with the problem of the second order 
aberrations of combinations of surfaces and systems, a problem which, so far as we 
know, has never been attacked from any general standpoint. Koenig and Von Eohr 
(Von Rohr, 'Theorie der Optischen Instrumente/ Cap. V.) give a development of a 
formula for the coefficients of first order and second order in the longitudinal spherical 
aberration, based on Abbe's method of Invariants, but so far as can be seen, no definite 
results are obtained for the second order terms. 

Dennis Taylor ('System of Applied Optics/ p. 67) gives a formula for the 
spherical aberration, developed in powers of the intercept made by the ray on the 
first principal plane, which includes terms of second order. But his formula, a 
particular case of those dealt with in the present paper, is limited to the thin lens, 
and no attempt seems to be made at anything like a general treatment of such 
aberrations. 

Another important object of the method to be described is to express the 



SPHERICAL ABERRATION FOR A SYMMETRICAL OPTICAL SYSTEM. 31 

aberrations in such a form that, in a combination of surfaces and lenses, the effect of 
a given surface or lens on the final result can be readily traced. This is fundamental 
for the designer, who usually proceeds to sketch out his system by Gaussian methods 
only, being guided therein by considerations of magnification, illumination, and field 
of view ; and then goes on to eliminate the resulting image defects, so far as he can, 
by bending the lenses, i.e., by altering their mean curvature without changing the 
focal length. In doing this he usually corrects one defect at a time, with the frequent 
result that, when, having corrected one defect by means of one lens, he proceeds to 
correct a second defect, he thereby causes the reappearance of the first. 

If the effects of any given lens, however, are made apparent in the final formula, it 
becomes a more manageable problem to devise variations which will keep any one 
defect invariant whilst others are being dealt with. 

§ 2. Notation. 

There is no general agreement among mathematical writers as to the notation 
employed in dealing with optical problems, and it will be convenient to state here 
the symbols we have adopted. They are a modification of a system due to 
Steinheil. 

The successive media, proceeding in the direction of travel of the light (from left to 
right in our figures), are denoted by even suffixes 0, 2, 4, &c, and the same suffixes 
affect the rays in these media, their inclinations, a , a 2 , a 4 , &c, to the axis, and their 
intersections I , I 2 , I 4 , &c, with that axis. 

The successive geometrical images will be denoted by the letter J, thus J , J 2 , J 4 , &c. 
The successive surfaces of separation will be denoted by the odd suffixes 1, 3, 5, &c, 
and the same suffixes will affect the centres of curvature, the intersections of rays 
with the surfaces, and the points where the axis crosses the surfaces. The latter 
will be denoted by the letter A and the centres of curvature by the letter C. 

Fig. 1 illustrates the use of this notation for two refracting surfaces. 

The radii of curvature are r u r 3 , r 5 , &c, and are to be considered positive when 
A 2n+1 G 2n+l is measured from left to right. 

The perpendicular from a centre of curvature on a ray is denoted by p and is 
affected by a double suffix, the first belonging to the centre of the curvature and the 
second to the ray. Thus p 12 is the perpendicular from the centre of curvature C x of 
the first reflecting surface upon the ray in the second medium.' Where there is no 
ambiguity the first suffix will usually be omitted. 

The refractive index will be denoted by n and affected by the suffix of its 
medium. 

Transverse magnifications will be denoted by M. The magnification produced by 
surface 1 will be denoted, as convenient, by M. 1 or M 02 ; by surfaces 1, 3 combined 
either by M l8 or M 04 ; by surfaces 1, 3, 5, combined either by M 135 , or M 0( ., and so on. 

F 2 
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The advantage of the double even suffix notation in this case is that we have a symbol, 
M 40 , for the magnification when light passes backwards through the system, the order 
of the suffixes being material. Where odd suffixes are used, we have to use l/M l5 
l/M 13 , &c., for the reversed magnifications. 

Ray magnifications will be denoted by M. These are the limit of the sine-ratio for 
small inclinations, thus M x = M 02 = L sin a /sin a 3 . M = M when the initial and final 

media are the same. 

With regard to inclinations, they will be treated as positive when the rays converge 
to the axis, as in fig. 1. The inclinations of the rays calculated by Gauss' process 
will be denoted by ft. Thus ft = a , tan ft 2 = tan a /M 02 , tan ft 4 = tan a 4 /M 04 , &c. 




/ <* P ^- 



o *o 



Fig. 1. 

We may also use angles y, calculated from a constant sine ratio, viz., y = a 0i 
sin y 2 = sin a /M 025 sin y 4 = sin a /M 04 , &c. 

Throughout much of the work we shall use the same trigonometrical function 
(tangent or sine) of the angles a. If the tangent is used, we shall employ the following 
abbreviations : — • 

$2n == ^ an a 2n ^2n = tan p 2n . 

If the sine is used, the meaning of g, t will be as follows : — 



q 2n = sin a 2n t 2n = sin y 2n . 

It will be found that many formulae remain unaltered, whichever of the two inter- 
pretations for q and t is used. 
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All distances parallel to the axis will be denoted by x (the attribution of the 
symbol being indicated in each case) and will be measured positively from left to 
right. 

The longitudinal aberration will be denoted by Ax, with the suffix of the medium, 
and will be reckoned positive when the point of intersection of the ray with the axis 
is to the right of the geometrical image. Thus Ax 2 = J 2 I 2 . This is opposite to the 
usual convention which is based on the fact that for positive or convergent lenses, I 2 is 
generally to the left of J 2 , ; but, in the first place, this is not universally true, and, in 
the second place, the convention adopted by us was found more convenient in handling 
the algebra. 

It is to be noted that, with the notation used, the well-known formula for a lens 

111 111. 

-* + - = — becomes — - - = — the distances u and v being measured in the same 

u v j v u j 

direction. 

The distances between successive refracting surfaces we denote bye, with the suffix 
of the medium. 

In dealing with a system, especially where the initial and final media are not the 
same, it is very convenient to use an " equivalent " Gaussian system, in which lengths 
parallel to the axis are measured in each medium in terms of a unit proportional to its 
absolute refractive index. 

If we denote the corresponding points in the equivalent Gaussian system by accents, 
we find that 

A / V _ A.xJ Q a / j/ __ A 1 J 2 a / a / ___ / __ C 2 

■**■ 1° — 3 ■**• 1° 2 — — ? JA. xJA. - s — C 2 — 5 

A' 3 j' 2 = A*L» , Ay, = ^J* &c . 

If then we denote the quantities — p , - A - — - , &c. , by -r , -j . . . ; fi,f 3 . . • may be called 

the focal lengths of the successive refracting surfaces. The equations connecting 
image and object in the equivalent Gaussian system are 



A iJ 2 A j J J x 

which is of the same form as the equation connecting image and object for a thin lens 
at AV 

Now bearing in mind that A 3 J 2 = AjJg— AiA 3 and therefore A^J^ = A^J^— -A! X A!^ 
it is easy to show that the effects of the successive refracting surfaces in the actual 
system can be obtained by compounding a corresponding set of thin lenses in the 
equivalent Gaussian system. By dealing with the latter, we get rid of the 
asymmetry introduced by the difference of initial and final index. Of course this 
applies only to the calculation of the geometrical images. We note that in the 
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equivalent Gaussian system the ray and transverse magnifications are identical and 
agree with the transverse magnifications in the actual system. 

Finally the intercept of the incident ray on the leading principal plane of a system 
will usually be denoted by y. This is taken by some authors as the argument of the 
development in series, but differs only by a factor from tan a or tan /3 2n . 

In many cases it will be convenient, in order to avoid unnecessarily large 
suffixes, to condense a system of surfaces or lenses, affecting quantities referring to 
the system itself with suffix 1, and the initial and final media with suffixes and 2, 
the paths in the intermediate media not being explicitly considered. 

§3. Singularities and Convergency. 

Consider any symmetrical optical system, of which PL and QM (fig. 2) are the 
initial and final refracting surfaces. Let F be the front focus of the system and 
UF V the caustic for backward- travelling rays which are parallel in the final medium. 



" — g.__ ■ » fc » 








Fig. 2. 

This caustic, as is well known, will usually be of the type shown in fig. 2, 
approximating to a semi-cubical parabola with a cusp at F , and, to fix ideas, we shall 
suppose the point of the cusp to be turned to the left. In the opposite case, an 
obvious modification of the argument will be found to lead to similar conclusions. 

Any ray in the initial medium, which touches this caustic, must emerge parallel to 
the axis after passing through the system. 

Let I be an object point on the axis behind F and sufficiently near to it for a real 
tangent to be drawn from I to the caustic and yet go through the system. This 
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will involve for I a positive ray magnification M exceeding some definite finite limit. 
The Gaussian image point J 2 will then lie a finite distance in front of F . 

Consider first a nearly paraxial incident ray I P. Such a ray will be refracted 
approximately according to the Gaussian law and will emerge at an inclination a 2 , 
where a 2 is nearly equal to a /M and both a and M being finite and positive, a 2 is 
also finite and positive. The ray emerges as QR, passing through a point I 2 , finitely 
different from J 2 . 

As a increases, a 2 at first increases with it, but as a reaches the value X, 
corresponding to the inclination of the ray I L which touches the front focus caustic, 
a 2 is again zero. Hence between those two values a 2 has at least one maximum, and 
for a given value of a 2 there are at least two values of a . 

Thus, within the range of values which are of practical importance, a is a many- 
valued function of a 2 having one or more branch-points, of which the one of least 
modulus corresponds to the first maximum of a 2 . 

Now, within the same range of values, all the aberrations must be given as single 
valued functions of a , since clearly there can only be one physical emergent ray, 
corresponding to one given physical incident ray. This statement, as we shall see, 
needs to be qualified when we are dealing with purely geometrical rays, but this need 
not affect the present stage of the discussion. 

In consequence, if any aberration be expressed in terms of a 2 — or of any trigono- 
metrical function of a 2 — that aberration must, in general, be a many -valued function 
of a 2 , having for its branch-point of least modulus the first maximum value of a 2 
mentioned above. It follows by a well-known result in theory of functions, that 
no Taylor's series in a 2 , or in sin a 2 , or tan a 2 , can be valid for values of a 2 
exceeding this modulus numerically. For such values the series will be definitely 
divergent. 

It is interesting to consider what happens when I is on the other side of F , so 
that we are dealing with a large negative magnification. In this case no real 
tangent can be drawn from I to the front focus caustic and the value of a , for 
which a 2 =0, is a pure imaginary. But here again, although we are now dealing 
with imaginary values, we get two values of a for a given (pure imaginary) value 
of a 2 , and, although no such maximum of a occurs in the purely real values, the 
modulus of the imaginary branch-point limits the validity of Taylor's series in a 2 as 
before. 

Thus there exists always a certain range, extending a finite distance (depending 
on the nature of the optical system) on either side of the front focus, within which 
no development of any aberration in powers of a 2 or of its trigonometrical functions 
(or, indeed, by similar reasoning, of any inclination of the ray, except in the 
original medium) is valid for the whole pencil of rays which actually traverse the 
system. 

Indeed, as the object point I approaches the front focus, it is clear that both X and 
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the maximum a 2 tend to zero, so that only an infinitesimal portion of the rays can be 
dealt with by the method of successive aberrations, i.e., by the Taylor's series. 

That the range of failure is by no means an unimportant one is shown by an 
example given by the authors in a paper read before the Optical Society in December, 
1918. In this example the system considered is a positive lens of unit focal length 
and thickness ^r, meniscus shaped, with curvatures 1 and 2*36, and its convex side 
towards the incoming light. For such a lens and magnification as low as 2, the 
critical value of a 4 is found to be about 4° 40', corresponding to a value of a of 13°, 
whilst the greatest practical value of a is 26°, so that in this case only about \ of the 
light going through the lens could be dealt with by series in terms of the emergent 
angle. From M = 2 to M = oo the conditions are still worse. 

As a matter of fact, it appears that in this case the range of magnifications, within 
which development in terms of the emergent inclinations is possible for all rays 
travelling through the lens, is restricted to a range lying somewhere between M = — 1 
and M = 1*5. This makes it clear that we cannot depend, in the calculation of the 
aberrations of an optical system, upon any series with the emergent inclination as 
argument. This is important, because from other considerations it would have been 
valuable to have been able to express the equation of the emergent ray in the 
form 

y + qx=f(q) 

where q is the inclination of the emergent ray, and to proceed to obtain successive 
approximations to the caustic by developing f (q) in powers. It now appears that 
this is not, in general, legitimate. 

We now come to the consideration of series proceeding by powers of a , or of its 
trigonometrical functions. Here the question of many-valuedness will not occur, 
except as follows. 

If we consider a ray impinging upon a spherical refracting surface, this ray, if 
produced, will meet the surface at a second point. Treating the problem from the 
purely analytical standpoint, this second point is also one at which refraction takes 
place, and thus, for the same a , there will, in general, be two values of a 2 , four of a 4 , 
and so on. a 2n will therefore, in general, be a multiple- valued function of a , and 
the aberrations will also be multiple-valued functions, and the branch-points of these 
multiple-valued functions will, as before, limit the convergency of the Taylor 
series. 

Now clearly two branches coincide whenever there occurs a grazing incidence ; 
and, therefore, if the system be so arranged (as it almost necessarily is) so that no 
grazing incidence is reached, there will be no real branch-points within the range of 
practical values. But this does not mean that the Taylor's series will necessarily 
be valid, for there might be imaginary branch-points. A very simple example will 
show how such branch-points can occur. 
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If I be a source of light placed in front of a plate of thickness c 2 and refractive 
index n> the perpendicular from I on the plate being the axis of the system, it is 
easily verified that the longitudinal spherical aberration 



J4I4 — C 2 

where sin a Q = n sin a 2 , so that 





1 __ tan a 2 
n tan a 0/ 

y 1— sin 2 a 

v n 

The branch-points here correspond to a = \ir or a 2 = j^ir, i.e., to grazing incidence at 
the first or second surface respectively. 

Clearly if n > 1, then, since sin 2 a rS 1, the second grazing incidence can never occur 
for real values of a . 

But if we take as our argument t — tan a , which removes the first branch-point 
to infinity, we find 



mammmmmm 





and this has imaginary branch-points where £ == + —77—5- — t . The radius of con- 

vergence of the Taylor's series in t is therefore given by t — /( 2 — -r, a value which 

does not correspond to any physical limitation of the rays. This applies to both the 
longitudinal and the transverse spherical aberrations in this case. 

The above example also brings out another important point ; for if in it sin a is 
taken as the argument, the branch-points are ±1, ±n; both of which correspond to 
definite physical limitations, viz., grazing incidence and total internal reflection, so 
that in this case the limitations of the Taylor's series are also the limitations of the 
problem. 

We see then that the validity even of the expansion in a may be limited by the 
existence of branch-points, and that the choice of the particular trigonometrical 
function in which we expand may exercise a considerable influence on the result. 

The limitation of the a developments due to branch-points will not, however, as in 
the case of the a 2 developments, lead to vanishing radii of convergence. There is 
always a finite region within which these developments may be used. In what 
follows, therefore, we have exclusively used a as argument. 

In dealing with the longitudinal Spherical aberration another limitation presents 
itself. We have seen that if a = X (fig. 2), a 2 = 0. It follows that the intersection 
of the emergent ray with the axis is then at infinity, or the longitudinal aberration is 
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infinite and afterwards changes sign. Thus the values a = ± X correspond to poles of 
the longitudinal aberration. These poles, being the singularities of least modulus, 
govern the convergence of the Taylor's series in this neighbourhood, and this radius 
of convergence tends to zero as the object approaches the front focus. It was a con- 
sideration of this difficulty which primarily led us to put the longitudinal aberration 
in a new form. 

This difficulty does not arise with the transverse spherical aberration. The poles of 
the longitudinal spherical aberration are due to the zeros of a , and on multiplying 

by — ^p— ° to get the transverse aberration, these poles disappear. 



M 



4. Summary of Method and Results, 



The general principle of the method employed was suggested by an attempt to fit 
an empirical formula to the longitudinal aberration of a lens for a certain range of 
curvatures and object and image positions. This empirical formula was discussed by 
the authors in a paper recently read before the Optical Society^ and was found to 
give, on the whole, a singularly good fit. Briefly stated, the formula is of the 
following type : — 

Ax== T^W (1) 

where t is the slope of the emergent " Gaussian " ray, so that t = t /M . A is the 
(known) theoretical constant of the first-order aberration, which is a quartic in the 
magnification, and B is a cubic in the magnification, the coefficients in which are 
determined empirically. This formula was found to give a good approximation, even 
when the magnification was high and we were working well outside the limits of con- 
vergency of the Taylor's series for Ax. 

If we consider any given object point, the longitudinal spherical aberration will be 
a function of t 09 that is, of t. Denoting it hyf(t) 9 the reasoning of the preceding section 
shows that f{t) is always one- valued for a finite (and generally quite considerable) 
range of t 9 but it is not regular, having poles at t = ± t, where r = tan x/M and 
becomes rapidly small as the magnification increases numerically. 

If, however, we write (l—t 2 /T 2 )f(t) = <j> (t) 9 <j> (t) is now limited only by the original 
branch-points off(t) and will, in general, have an adequate radius of convergence. 
We may therefore expand it in a Taylor's series, and we get for/ (t) the form 

A ___ > /,\ __ at + bt + ct 4* . . . /q\ 

J \ I l-~-t 2 /r 2 V J 

* Baker and Filon, " On an Empirical Formula for the Longitudinal Spherical Aberrations in a Thick 
Lens," < Proceedings of the Optical Society/ December, 1918. 
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If the series in the numerator converges rapidly, it will be sufficient, provided a is not 
near zero, to stop at the first term, and we get as an approximate formula 

A *=T^P <»> 

which is of the same form as (l). 

We have necessarily a = A, and if the two formulas are to tally we should have in 
addition B = — 1/t 2 . 

The formula in the form (3), however, is not rigorously correct to the second order of 
aberrations inclusive, unless b happens to be small. If we wish to retain second order 
terms complete, we have to use 

a — atf2 + ^ (a\ 

~~ fZ~fO v 4 / 

and this can be written, to the same order of algebraic approximation, in the form 

Ax = at 2 /{l-{l/r 2 + b/a)t 2 }, . (5) 

provided again a is not zero. 

If this form (5) is adopted, then the B of the empirical formula should be \\T+b\a. 
But if this is done, the formula suffers from two defects ; (i) it fails whenever a is 
near zero ; (ii) it does not give exact compensation for the poles in the critical range 
for M large. 

The further question then arose : how far are formulae of type (4) or (5) suitable 
for dealing with combinations of surfaces or lenses ? An important guiding considera- 
tion, in all work of this kind, must be the relative simplicity of the formulae in 
passing from a single surface or lens to a combination, and whether these formulae are 
suitable for tracing the effect of individual surfaces or lenses upon the final result. 

We have ultimately been led to the conclusion that no single formula can satisfy 
completely the three ideal requirements, viz. : (i) exact agreement with development 
as far as the second order inclusive ; (ii) simplicity in dealing with combinations ; 
(iii) exact compensation of the poles in the critical range of M. 

The method finally adopted satisfies conditions (i) and (ii). It only satisfies (iii) 
approximately. Numerical calculations show that numerically the approximation is 
adequate in the case of a lens or a simple surface. In the case of more complicated 
systems we have, as yet, no numerical data. 

The first part of the investigation deals with the single refraction. It is there shown 
that the longitudinal aberration can be put into the form (l), i.e., Ax = At 2 /(l + Bt 2 ), 
where the formula is correct to the second order inclusive. 

We also find, for the inclination of the emergent ray, the formula 

q = t(l + Bt 2 )l(l + Ct 2 ), (6) 

U 2 
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which is correct to the first order when q and t are tangents and to the second order 
when q and t are sines. * 

In the above A, B, C are polynomials in M of degrees 4, 3, 2 respectively, so that 
the empirical formula is well justified for the simple refracting surface. In this case, 
too, it is possible to calculate t directly, and, in fact, a simple geometrical construction 
is given for it. When this is followed for varying image-positions, it is found that 
outside a certain range of M, the r so obtained becomes irrelevant, and that, in fact, 
if the correct factor 1— t 2 /r 2 is retained in the denominator of Ax, although it 
improves the fit by removing singularities in the range round M = oo 9 it introduces 
entirely fictitious singularities in other and important parts of the range, and makes 
the formula worthless. 

A good deal of light is thrown upon the problem when it is found that, if we 

1 
develop — in descending powers of M in the neighbourhood of M = oo ? the two 

T 

leading terms are discovered to be identical with the two leading terms of the 
cubic B, previously obtained. This makes our B approximate more and more closely 

to -^precisely as the effect of the denominator term becomes more important, and it 

T 

is this fact which is the key to the numerical value of the method. 

We then proceed to show how the constants for a combination of the two systems 
can be obtained from the corresponding constants of the individual systems. In 
doing this it appears that, so soon as we pass from the single refracting surface to 
the lens, a new constant is introduced into the formula, which now takes the form 

■ *-w (7) 

where A and B are of the same form as before, but E is now a polynomial of 

degree 6 in M. In the case of a lens the term E£ 4 is found to be, in general, of 

small importance, which accounts for the good fit of the empirical formula. 

The formula (7) for a combination holds good to the second order inclusive, and B 

1 
agrees with — , when M is large, as far as the leading term only. For numerical 

purposes, however, a correction is discussed, which is very readily applied, and which 

makes the two leading terms in B agree with the two leading terms in -§, as in the 

T 

case of the single refracting surfaces. 

The formulae for combining two systems take comparatively simple forms ; the A, 
B and C for the combination are expressed as linear functions of the A's, B's and C's 
of the components, and the E as a lineo-linear function of the A's, B's and C's of the 
components, each term involving a product of which one factor belongs to one 
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component of the combination and the other factor to the second. In addition the E 
for the combination involves linear terms in the E's of the components. 

These results are found to hold good in the more general case of the combination 
of three or more systems. It will follow that if the constants A, B, C, E are 
tabulated for lenses of various curvatures, the effect of the combination can be traced 
relatively easily and the aberrations corrected, so far as possible, by suitably bending 
the lenses, while keeping the general arrangement and the magnifications the same. 

Explicit values of the constants for the single refracting surface and a single thick 
or thin lens have been obtained and are tabulated for reference, so as to be available 
for eventual computation of the required tables. We have also given some numerical 
values for a single lens, and a numerical test of the accuracy in this case, which works 
out at about ^q of the total aberration for the range of cases taken. 

The corresponding formulae with sin y instead of tan /3 as argument are discussed, 
and it is shown that the equations of combination are of the same form as before. 

Certain invariant relations between the coefficients in A, B, C, E are developed, 
which enable various calculations to be simplified and in particular to determine these 
constants for a system reversed, when they are known for the direct system. This 
will generally halve the work of tabulation. 

§ 5. The Single Refracting Surface, 

Using the general notation described in § 2, consider refractions at a single refracting 
surface. 

Let \^ , \fs 2 denote the angles of incidence and refraction, so that \fs = CJPJo, \[s 2 = 
CiPJg (fig. 1). 

Let Gilo = X == x Q , CJg = X 2 = x 2 + Ax 2 , CjJg = x 2 we then have the set of refraction 
equations 

sin ^ = p /r 1 = x sin ajr 1 (8) 

sin y\r 2 = j9 2 / r i = X 2 sin a 2 /r 1 (9) 

2 A' 9 " — ^0 /^0 * * * * * * * * * * «IX\_// 

a 2 — ot, = V r o~~"V r 2 . .(ll) 

Let £ = -A/J'oj ^ 2 = A/iJ ; 2 in the " equivalent" Gaussian system (see § 2). Then 

n £o = Xo + r l9 n 2 £ 2 = x 2 + r u (12) 

and we find, using the first approximation when a , &c, are small 

!/&-!/&= Wi- • • • • • ( 13 ) 
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If Mj = La /a 2 , which we shall call the ray magnification, we have 

a ->0 

The transverse magnification M^ is given by 

Mi = ^o = ^ = ^/% • ....... (15) 

We can also express (13) in another well-known form, namely, 

n /x 2 -n 2 /x = l/f u (16) 

whence, using (14), we obtain 

x = n 2 f 1 (l-M 1 )/m i ^ 

^ = ^0/1(1-^) J ' 
Again, from (9), (8) and (10) 

X 2 = £> 2 /sin a 2 = 7i £C sin a /n 2 sin a 2 = w j^ (l — M x ) sin ajMj sin a 2 

= # 2 sin a plL 1 sin a 2 , (18) 

and the longitudinal aberration 

Ax 2 — X 2 — x 2 = x 2 (sin oLofMLi sin a 2 — l). ...... (19) 

The corresponding longitudinal aberration in the equivalent Gaussian system is 

found from 

A£ 2 = Ax 2 /n 2 — (w 0t /i/w 3 ) (l — Mj) (sin ao/M^ sin a 2 — l) (20) 

Now from (ll) 

sin a 2 = sin a cos i/r cos i/r 2 + sin \[s cos a cos i/^ 2 — sin \fr 2 cos a cos ^ 
+ sin a sin i/r sin i/r 2 , 

whence, using (8), (9) and (10), 

sin a 2 /sin a = { l—(x /r l ) 2 sin 2 a p { 1 — (ngxjn^) 2 sin 2 a }* 

+ (^oM) (J -sin 2 a P { l-in.xjn.r,) 2 sin 2 a p 
-(WVi) { 1 -sin 2 a }* { 1 -(xo/n) 2 sin 2 a }* 
+ (v 2 /n 2 n 2 ) sin 2 a , 

and developing this in ascending powers of sin a , we obtain, retaining only terms of 
fourth degree 

sin a 2 /sin a = 1 -f cc ( n 2 —n )/n 2 r 1 —^ : P sin 2 a 

~|P sin 4 a {(^ 2 + ^ )X 2 + ^ 2 (r x -x ) (n^ + n^)} ln 2 2 r x 2 , 
where 

P = (l-njnv) (»oM) (1 +^ /n) (l~v /vi) ? 
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and is the quantity whose vanishing gives the aplanatic points and must therefore be 
a factor of every coefficient after the first in the development of sin a 2 in powers of 
sin a . 

If we write for shortness 

Q = (l + n /n 2 ) 2 x 2 /r 1 2 + (l ~ x o/n) (1 +VoMn), 

and remember that 

1 + % (w 2 -ti )/%n '=' 1/M 1? 

we find 

sin a 2 /sin a = l/M^— |\P sin 2 a /(l — |-Q sin 2 a ) 

= MrMl-ftPMi + iQ) sin 2 a }/(l-iQ sin 2 « ) 

= W i - 1 {l+Bam 2 ot /M 1 2 }/{l + C8iii 2 a l) /JIL 1 2 } ... (21) 

correct as far as the second order inclusive, where 

B = -IPM^-IQM, 2 , C = -iQM, 2 , 

from which, after some reductions 

C = -i(n 2 -n )- 2 {« + n ra 2 + O-3 (n 2 + n 2 2 ) M, + 3 {ni-n^+nfiW. 2 }, (22) 
and 

= i (n 2 -n )- 2 { -(n 2 2 -n n 2 +n 2 ) + (n 2 -n fM 1 -(n 2 2 +n 2 -5n n 2 ) M 1 2 -2w w 2 M 1 3 }.(23) 
Keturning to equation (19) and using (21) 



where 



Ax 2 = x 2 (C -B) Mr 2 sin 2 « /{ 1 + BMr 2 sin 2 « } 

= w 2 / 1 AMr 2 sin 2 a /{l + BMr 2 sin 2 a }, (24) 

A = x a (C-B)/r h f 1 

= -I (»2-w )- 2 (»oK) (1 -Mi) 2 (nj-noM,) (« -n 2 M x ) 

w \ 2/ / x r-n n 2 + (n 2 +n ) 2 M 1 -2(n 2 2 +n n 2 +n 2 )M. 2 l 

— i\n 2 — n )~ (n /n 2 )< , s ' >. (25) 

All the above formulae are correct to the second order of aberrations inclusive. We 
note that A, B and C are polynomials of degree 4, 3 and 2 in the magnification 
respectively. 

If we express the aberration in terms of tangents instead of sines we have at 
once 

Ax 2 = w^AMr 2 tan 2 a 3 /{ 1 + (B + M/) Mr 2 tan 2 « } 

= w 2 / 1 Atan 2 /3 2 /{l+Btan 2 /3 2 } . (26) 
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where 

b = B+n* 

= i(n 2 ~-n )- 2 {~ (n/~n n 2 +n 2 ) + (^ 2 ~^ 

When, however, we come to develop a formula for tangents, similar to (21) for 
sines, it is found that, in the series for tan a 2 /tan a , we do not have all the coefficients 
after the first vanishing together ; for, even at the aplanatic points, the tangent ratio 
is not constant. 

We can, indeed, write 

tail a^tana = M 1 - 1 (l + BMr a tan 3 Oo)/(l + CMr 2 tan 2 oo), . . . (28) 

and choose the coefficients B and so that the developments shall agree as far as terms 
in tan 4 a inclusive. This can, in general, be done in one way only. But we then find 
that B and are no longer integral functions of the magnification. Their infinities 
have to be taken into account, and generally the method becomes complicated and 
unsatisfactory. 

From other considerations, however, it appears that since the zeroes of a 2 must be the 
same as the poles of Ax 2 , the B in (28) must be the same as the B in (26), and this 
will fix G as follows : — 

tan 2 a 2 = sin 2 a 2 f(l —sin 2 a 2 ) 

M^ 2 sin 2 a { 1 +B sin 2 ajM^} 2 



{ 1 + C sin 2 a /M 1 2 } 2 ^Mr 2 sin 2 a { 1 + B sin 2 ajW^ 2 } 2 . 

= Mr 2 tan 2 a { 1 + (B ±W) Mr 2 tan 2 * } 2 

(l+tan 2 a ){l + (C + M 1 2 )Mr 2 tan 2 a } 2 ~Mr 2 tan 2 a {l + (B + M/)M 

substituting from (21). 

Hence, taking the square root and developing the denominator in powers of tan /3 29 
i.e., of M^"" 1 tan a , 

, n o __ ( 1 4- B tan 2 ft 2 ) 

21 ^ ~ l+tan 2 /3 2 (C+fM/-|^ 

= (l + Btan'^^l+Ota^^ + Dtan 4 ^), (29) 

where 

C - C-f fM/-! 

D = C(l+M 1 2 )+|M 1 4 + fM 1 2 -i~2B L. . . . . (30) 

= C(l+M 1 2 ) + J(l- M i 2 )(l + 3M 1 2 )-2B^ 

(29) is now correct as far as the second order inclusive. 

If we only require the tangent ratio correct up to the first order of aberrations, 
we have the formula 

tan a 2 = tan /3 2 ( 1 + B tan 2 /3 2 )/( 1 + C tan 2 /3 2 ) ..... (31) 
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The value of C, when written out fully, is given by 

C = f (n 2 -n Q )~ 2 {-{n 2 2 -n n 2 + n 2 ) + [n* + n 2 2 ) M 1 + (w 3 -3w n a + n 2 2 ) HL 2 }. . (32) 

§ 6. 7%6 Convergency Factor and the Singular Inclination for a Single 

Refracting Surface. 

Having now obtained expressions (24 and 26) for the longitudinal spherical 
aberration, which are correct to the second order of aberrations, when expansion 
in powers of sin a or tan a is legitimate and rapidly convergent, we have now to 
enquire how far the same expression remains valid as M increases, in which case we 
know that B or B increases without limit and the convergency fails, even for 
comparatively small values of a . 

Here it will be convenient to introduce two definitions :-*- 

I. We shall call singular inclination the value X (see § 3) of a for which the 
emergent ray is parallel to the axis. 

II. The factor 1 — sin 2 a /sin 2 X (or 1 — tan 2 a /tan 2 X if we are dealing with tangents) 
we shall call the convergency factor. If we multiply Ax by the convergency factor 
we remove those singularities of Ax which are instrumental in causing critical failure 
of convergency. 

To find the singular inclination and convergency factor for a single refracting 
surface, we have to find when a 2 = 0. 

Going back to the fundamental equations (8) to (ll) we have a 3 = when a = X, 
where 

X = ^2-^0 

which leads to 

sin X = sin \fr 2 cos \fs — sin yfr cos \Js 2 

= (n x sin X/n 2 r x ) </{ 1 - (x sin X/r,) 2 } - {x sin x/r,) ^/{ 1 - {n x sin X/n 2 7\) 2 }. 

Hence, either sin X = 0, which obviously refers to the axial ray, a trivial and (for our 
purpose) irrelevant solution, or 

Tl /x = {n /n 2 ) \/ {I -{xo sin Xfa) 2 }- ^{l-in^sin Xjn 2 r x ) 2 }. . . . (33) 
On rationalising (33) leads to 

4n 2 sin 2 x/n 2 2 = ±r x 2 jx 2 - ( 1 + rflx* - n^jn, 2 ) 2 

= - ( 1 + rfx, - njnj) ( 1 4- r fx + n /n 2 ) ( 1 - r x jx + n /n 2 ) ( 1 - r,/x () - njn 2 ) .(34) 

This gives the singular inclination. 
If we write 

K = (1 +r,/ir -?i /?* 3 ) (1 +r,/aj + ?io/w a ) ( L -f,/o: -fnj%) {l-r } /xo-nJri 2 ) f . (35) 

VOL. CCXXI. — A. H 
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it follows that 



1 4- 4n 2 sin 2 a /n/R . 



(36) 



is the required convergency factor. 

If the formula of §5 are to get accurately over the failure of convergency, this 
convergency factor should be identical with 



that is, we should have 



1 + B sin 2 ao/Mj 2 , 
B = ±n 2 M 1 2 /n 2 2 'R, ; . (37) 



which, when written out, becomes 

B = tooVI/ {l-^iY/(n 2 -n ) 2 (1-2M!) (n 2 +n -2nJSL) {n^n.-Zn^) . (38) 

This does not agree with the previously found value for B, being of fractional form 
in M 1# It does lead to B becoming infinite of the order M^ 3 when M 1 tends to infinity, but 
it indicates an infinity of B (and therefore a critical failure of convergency) at three 
other places, namely when M x = |y|- {n -\-n 2 )/n , \ (n + n 2 )/n 2 , at none of which does a 
failure of convergency really occur, as can readily be verified. 

The reason for this is made clearer by geometrical reasoning as follows : — 

Let PiQi (fig. 3) be a ray which is parallel to the axis in medium 2. To make the 
figure easier and the quantities dealt with positive, the refraction has been taken 
from a denser to a rarer medium, so that n /n 2 '>l* 




Fig. 3. 

In the triangle IoC^ of fig. 3 we have sin \Jr 2 j sin \fs = loP^^. 

But 

sin \j/ 2 ] sin \Js Q = n /n 2 ; hence 1^ = n^\n 2 
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The point P x and singular inclination X can therefore be constructed geometrically as 
follows. 

With I as centre and radius ^.CJo/^ describe a circle meeting the refracting 
surface at P x ; CJqPi is the angle X required. 

This angle X approaches zero, that is, we get a critical failure of convergency, when 
the two circles approach contact at A^ The limiting case is, therefore, when 
AJo = n . CiIo/^2 or I divides A^ externally in the ratio n :n 2 . 

When this happens r 1 + x = n x /n 2 leading to M x = oo ? a case of true critical 
failure. But clearly, by symmetry, we get a precisely similar result when P1Q1 is 
due to a ray entering the surface at Q l9 and travelling backward through medium 2. 
In this case the limiting position of I divides B^ externally in the ratio n Q : n 2 and 
is defined by ri—x = —n x /n 2 or M^ = |-. 

This case would correspond, analytically, to a 2 = 7r, and the corresponding equa- 
tion (33) would become 

-n/a^ = (n Q /n 2 ) *J{ 1 - (x sin X/n) 2 } - a/{ ! ~-(n x sin X/^n) 2 }. 

Now if we examine (34) we find that in the process of clearing roots, n/^ appears 
squared in the final result, which accordingly includes both a 2 = and a 2 = it. If we 
write Tifx 2 = u 9 then we should really write equation (33) in the form 

y/u = (n /n 2 ) a/(1 — sin 2 X/w)— \/(l — n o 2 sin 2 \/n 2 2 u), . . . . (39) 

and the two cases are discriminated by assigning to y/u one or the other sign. One 
of these cases is necessarily irrelevant since refraction at the posterior surface of the 
sphere is physically excluded. 

Further, if we consider the other two values which make R = 0, viz., 

M x = (n + n 2 )/2nQ and M x = (n -\-n 2 )/2n 2 , 
they correspond to 

and 

Vi Xq = 71qXq/71 2> 

i.e., to positions of I in which it divides AA and Bid internally in the ratio n^':n 2 . 
But these belong geometrically to the limit of Cases in which the incident and 
refracted rays lie on opposite sides of the normal, i.e., to a negative refractive index. 
And indeed they are obtained from the two previous points by reversing the sign of 

n /n 2 . 

Here again, examination of (34) shows that n /n 2 appears squared in it. Therefore 

(34) includes the cases in question. These* however, may be obtained analytically by 

changing Yr , or ^ 2 , into its supplement, i.e., by reversing the sign of cos ^ or cos^ 2 

or by changing the determination of the sign of one or other of the square roots on 

the right-hand side of (39). 

H 2 
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The cases are discriminated by the vanishing of these square roots, which occurs 
when cos yjs or cos \fs 2 = 0. 

It appears, therefore, that the vanishing of the factors in the denominator of (38) is 
wholly irrelevant, and, if we adopted for B the value given on the right-hand side 
of that equation, we should thereby be introducing, in the neighbourhood of Mi = Jr, 
(#o + n 2 )/£n Qr {n + n 2 )/2n 2 entirely irrelevant singularities, which would make the 
formula worthless. 

The question arises, what is the range of values of M x for which the equation 

4% 2 sin 2 \/n 2 2 = — E 
is valid and legitimate ? 

If we start from M x = oo , which corresponds to a real case, the signs of the square 

roots in (39) are well determined, and the correspondence between sin 3 X and M 2 is 

unique and definite and can be continued until we reach a point where one case 

passes into another. These cases we have found to be the branch-points of the 

three square roots, namely :• — 

u = 0, cos \fs = and cos \/s 2 = 0. 

u = leads to x = oo or M x = . . . . . . . . (A) 

cos ^ = leads to sin 2 X = u 9 or, using the first form of (34) 

kn 2 u\n£ = 4^— (l+^~ n 2 jn 2 ) 2 , 

i.e., {l—u—n 2 /n 2 2 ) 2 = 0, that is u = l—n 2 /n 2 2 , leading to 

x Q = ±r x {l-n 2 /n 2 2 )S and M x = {l±(n 2 -n )/^{n 2 2 -n ')}- 1 ...(B) 

cos \{s 2 = leads to sin 2 X = nfujn 2 , that is, to 

u = nflnf-l, x = ±r x (n 2 /n 2 2 -l)~\ M x = {i±(n 2 -n )/ x /(n 2 -n 2 2 )}- 1 . . (C) 

If n > n 29 both values of M x given by (B) are imaginary. The values given by (0) 
are both positive, HL X = {l + (w a — Wo)/^^ 3 — n/)}- 1 being the greater. 

The range over which we can travel without ambiguity is, therefore, from 
Mi = -f go to M x = {I + t n 2-n )/ v / {n ( 2 -n 2 2 )}- 1 and from M x = - ■ oo to M, = 0. 

If n <n 2 , the values of M x given by (C) are imaginary, those given by (B) are 
positive, and M x — {l— (^ 2 ^%)A/(^/— V)}" 1 is the greater, so that the range of 
validity is from NL 1 = + oo to HL X = {l—(n 2 —n )/^(n 2 2 —n 2 )}"' 1 and from M 1 =, ~ oo to 
M x = 0. 

Within this range (T + 4w 3 sin 2 a /^ 2 2 R) is the correct convergency factor; outside 
this range it is irrelevant. 

It is clear, then, that we cannot find a single formula for the convergency factor, 
which will hold for all values of the magnification. 
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Further, if the factor ( 1 + in 2 sin 2 a fn 2 H) is introduced into the denominator 
of Ax 2 , we no longer obtain expressions of the simple type (24) and (26), and 
endless complications are introduced when we come to consider a compound 
system. 

Can we make our expression B given by (23) give a tolerable approximation to 
(4:n 2 M 1 2 /n 2 2 R) for those regions where the denominator factor is really needed, namely 
for M^ large, positively or negatively ? 

To get the answer to this question we develop (^n/M^/n/B,) in descending powers 
ofM x . 

This is found to be (the most rapid method is to break up first into partial 
fractions) 

i (%-^ )~ 2 f -2?i n 2 M 1 s -(n 2 -5n n 2 +n 2 2 ) M/ 

^ -i{n 2 -n o y + nzW^il2n^ V. (40) 

+ terms in l/M 1? &c. 

If we now compare (40) with (23) we find that the most important terms when M L 
is large, namely those in M/ and M^ agree in the two expressions. 

We may, therefore, take it that the approximations (21) and (24) which we have 
seen hold good to the second order when expansion in series is convergent, will 
probably not be numerically very far out when M x has a large value, in which case 
the normal method of development cannot be used. 

It is important, at this stage, and to justify the above assertion, to consider a few 
numerical examples. 

Tables I. and II. give the values of Ax 2 and sin a 2 for a single refracting surface, 
calculated for a number of values of M x and two inclinations in each case. The 
inclinations are fixed from the perpendicular distance tar of A x from the incident ray. 
This, for moderate inclinations, is sensibly the same as the intercept made by the 
incident ray on the principal plane, tar has been given the two values 0*5 and 0*25 
in every case, except for HL l = 2 where tar = 0*5 leads to a physically impossible value. 
In this case tar = 0*25 and tar = 0*125 have been used to define the ray. 

In each case four values have been computed (l) the correct one, from trigono- 
metrical calculation ; (2) the values given by formulae (21) and (24) — these are shown 
in the column headed '" fractional formula" ; (3) the values obtained by expansion in 
series, up to the optician's first order of aberrations inclusive, that is including sin 3 a 
in the development of Ax 2 and sina 2 — these are shown in the column headed •" first 
order" ; (4) the same series carried to the second order of aberrations inclusive, i.e., 
to the terms involving sin 5 a — these are shown in the column headed " second order." 
It should be noted that these first. and second order approximations are the most 
accurate that can be obtained, much more so than more usual ones, proceeding in 
powers of sin a 2 or tan a 2 . 
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Table I. — Values of Ax 2 for Single Refracting Surface. 



M. 


£7. 


First order. 


Second order. 


■ ; ■'■■ , ■ ■ 

Fractional 
formula. 


True. 


10 


0-5 


-33-35293 


-107*6252 


+ 27-18563 


+ 26-66794 


10 


0-25 


- 8-33823 


- 12-98025 


-18-81008 


-18-85273 


2 


0*25 


- 0-250000 


- 0-378906 


- 0-516129 


- 0-527526 


2 


0*125 


- 0-062500 


- 0-070557 


- 0-071749 


- 0-071781 


0-5 


0*5 


- 0-015625 


- 0-016541 


- 0-016598 


- 0-016611 


0-5 


0*25 


- 0-003906 


. - 0-003963 


- 0-003964 


- 0*003963 





0-5 


- 0-166667 


- 0-174769 


- 0-175183 


- 0-175809 





0*25 


- 0-041667 


- 0-042173 


- 0-042179 


- 0-042189 


- 1 


0-5 


- 1-000000 


- 0-947500 


- 0-950119 


- 0*956680 


- 1 


0-25 


- 0-250000 


- 0-246719 


- 0-246761 


- 0-246838 



Table II. — Values of Sin a 2 for Single Refracting Surface. 





73. 


First order. 


Second order. 


Fractional 
formula. 


True. 


10 


0-5 


+ 0*0250865 


+ 0-0454644 


+ 0-0576348 


+ 0-0610770 


10 


0-25 


-0 


•0078936 


-0-0072568 


-0 


•0071911 


-0-0071828 


2 


0-25 


-0 


•2187500 


-0-2065430 


-0 


•1987180 


-0-1978219 


2 


0-125 


-0 


•1210938 


-0-1207123 


-0 


•1206710 


-0*1206691 


0-5 


0-5 





•2539063 


0-2541962 





2542195 


0-2542230 


0-5 


0-25 





1254883 


0-1254974 


o- 


1254975 


0-1254974 





0-5 


0' 


1805556 


0-1817130 


0' 


1826667 


0-1827294 





0-25 


0- 


0850694 


0-0851267 


0- 


0851286 


0-0851291 ' 


- 1 


0*5 


0' 


1250000 


0-1299375 


o- 


1311526 


0-1314354 


.- 1 


0*25 


0-0531250 


0-0532793 


0-0532873 


0-0532884 



It appears from the above that the fractional formulae are not merely equal, but 
appreciably superior to the second order formulae, and this not merely in cases such 
as those of the three first entries in Table I. , in which the convergency of the series 
for Ax 2 is either absent or slow, but in every case where the fractional or second order 
formulae differ sensibly from the true value. (Clearly a divergence of 1 in the last 
place cannot be claimed as significant, for the last figure in Tables I. and II. is probably 
not correct within ±2, in some cases.) An estimate of the range of the formula can 
be obtained from the fact that in the cases, vs = 0*5, M x = 10 and 2, the angles of 
incidence were 52° 34' and 48° 35' respectively, and, for the other values, angles of 
incidence of 20° and 30° are quite common. , 

In view of this the accuracy of the results is surprising and, from the point of view 
of the further applications of the method, most encouraging. 
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§ 7. Combination of Two Systems. 

Call the systems 1 and 3, and the initial, intermediate and final media 0, 2, 4. 

fu fs are the focal lengths of the systems, as defined in § 2. M b M x are the trans- 
verse and ray magnifications in the first system, M 3 , M 3 in the second system. 

M 3 +AM 3? M3 + AM3 refer to the transverse and ray magnifications in the second 
system when I 2 , the true intersection of ray 2 with the axis, is taken as the object 
point for the second refraction (instead of J y , which refers to transverse and ray 
magnifications M 3 , M 3 ). 

Using the notation of §§ 2, 5, we assume 

Ax 2 = n 2 f 1 (A 1 t 2 2 +E 1 t 2 *)f(l+B 1 t 2 2 ) (41) 

q 2 =t 2 (l + B l t 2 2 )/{l + G 1 t 2 2 ) (42) 

where q 9 — , lo,}j 2 = , ^? k and B l5 Oi have suitable forms according as sines or 

* 2 tan 1 " J 2 [tan/3J 5 15 & 

tangents are considered. The constant E x is zero if the systems reduce to single 
refracting surfaces. Its form in the more general case will be discussed later. 

If we denote by A^ that part of Ax± which is due to Ax 2 and by A 3 a? 4 that which 
is introduced by the aberrations proper to the system 3, 

A 1 x 4c =-n,f 3 AM 3 

where AM 3 is obtained from Ax 2 by means of 



leading to 



Thus 



Ax, = n 2 f s {1/(M 8 + AM 8 )-l/M 8 }, 

AM 3 = -M 3 2 Ax 2 /(n 2 f s +M s Ax 2 ). 



A x x, = nJM* (A^ + E^/VIl +«2 3 .(B 1 + M a A 1 /i//8)} (43) 

Again 

. _ nj s {(A 3 + AA 3 ) ql (M 3 + AM 3 )- 2 + (E 3 + AE 3 ) q.f (M 3 + AM 3 )~'} 
^ 4 " l + (B 3 +AB 3 )g 2 2 (M 3 +AM 3 )- 2 

= » 4 / 3 [{A 3 M 3 - 2 + AM 3 d(A 3 M 3 - 2 )/dM 3 }g 2 2 +E 3 ^M 3 -*]/(l + B 3 ^M 3 - 2 ), (44) 

retaining only terms of second order in t 2 3 . Writing now for q 2 2 in the above its " first 
order" equivalent i 2 2 + 2 (Bj—Cj) f 2 4 , we have 

Axjrii = A 1 xjn i + A^cjfii 

{/ 1 M 3 2 A 1 +/ 3 A 3 /M 3 2 } tf+tf [/ 1 {B 3 A 1 M 3 2 M 3 - 2 +E 1 M 3 2 -A 1 M 3 2 d (A 3 M 3 ~ 2 )/<M 3 

.._ +A 1 A 3 M 3 M 3 - 2 }+/ 3 {B 1 A 3 M 3 - 2 +2A 3 (B 1 -C 1 )M 3 - 2 +E 3 M 3 ^} l . . 

{l+^^ + A^/^Hl + B^/M/} ••••. V40| 
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Remembering that t 2 = M 3 £ 4 and retaining only the first two terms of the denominator 
product 

Aa?A = (/iaAi^+MA 4 )/(l f W" (46) 

where 

/ 13 A 13 =/ 3 A 3+ /AM 3 2 M 3 2 .......... (47) 

B 13 =B3+B 1 M/+A 1 M 3 M 3 y i //;. ....... (48) 

/ 13 E 13 =/ 3 E 3 +/ 1 E 1 M 3 2 M 3 * 

+ /AM 3 2 {B 3 M 3 2 +A 3 M 3 -M 3 2 M3 2 d (A.^L :i ~ 2 )/dM s } 

+ / 3 A3M/{3B;-2C 1 } (49) 

. q 2 (M 3 + AM,)- 1 { 1 + (B 8 + AB 3 ) q* (M 3 + AM 3 )~ 2 } 



Again 



Q* 



l + (C 3 + AC 3 )g/(M 3 + AM 3 ) 



-2 



and retaining only terms of order t 2 3 

g 4 = q 2 (M 3 + AM3)- 1 (l + BA a M 3 - 2 )/(l + C 3 ^M 3 - 2 ) 

= ^m 3 - i (i-am 3 /m 3 )(i+b^)(i+b 3 ^m 3 - 2 )/{(i+g 1 ^)(i+Ca 2 m 3 - 2 )} 

= t 2 M 3 -\ ( 1 + B^ 2 + B 3 ^ 2 M 3 - 2 + tJAM^mftl + W + C 3 « 2 2 M 3 - 2 ) 

= t i (l + B ls t*)/(l+C ls 0, ...... f (50) 

where B 13 has the value given by (48) and 

The equations (47), (48), (49), (51), give the constants for the combined system in 
terms of those for the components. It may appear at first sight as if the choice of 
the constants B 13 and E 13 had been arbitrary, for clearly, if X be any quantity, 

{/lsA^z+aCe^+xa^^i/ii+^b^+x) */} 

will give a development equally valid to the second order. But, if we do this, and 
we wish to preserve the simple character of the relation (48) giving the B for the 
combination, X will have to be a linear function of A 1? A 3 , B 1? B 3 . XA 13 must then 
necessarily contain terms of one or other of the forms A^, A x 2 , A 3 B 3 , A 3 2 . Thus 
the new E will contain such terms and will no longer be of type (49) which is linear 
in the aberration coefficients of each system tak^n separately. Thus the lineo-linear 
type of equation for E 13 requires X = 0. 

We note also that the equations of combination are identical in form, whether we 
are dealing with the sine or the tangent of the inclination as argument. 
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§ 8. Nature of the Quantities A, B, C, E in the General Case of any System, 

In the case of the single refracting surface we found that A, B, C were polynomials 
of degrees 4, 3, 2 in M, and that E was identically zero. 

In addition, for such a surface, equations (23), (25) and (27) show that the 
coefficient of M t 4 in A is n /n 2 times the coefficient of M x 3 in B or B. This may be 
otherwise stated in the form : — 

A 1 -n M 1 BjM 2 , i-&> Ai-MJBi (I) 

is of the fourth degree only in appearance and reduces to an expression of the third 
degree in M 1 or M x . 

The same holds good for A 1 — MiBj, so this result is independent of whether the 
sine or tangent is taken as argument. The same remark applies to all the results of 
the present section and to the other invariant relations shortly to be proved. We 
may therefore conveniently state it here once for all. 

If we now refer to the equations (47), (48), and remember that in any combination — 

M 1 3//3 = M 3 // 1 3-1// 1 (52) 

and 

l// 1 M 13 =l// 13 M 1 -l// 3 (53) 

with the corresponding equations 

n ^is/f = n 2 M. s /f s -nJf (54) 

4i M i3 = 4i3Mi-43 • (55) 

and the obvious conditions 

M 13 = M 1 M 3 ; M 18 = M^, 

we note first that, if A 3 is a quartic in M 3 it is also a quartic in M 13 or M 13 , and that 
if A l is a quartic in M 1? it becomes, on multiplication by M 3 2 M 3 2 , a quartic in M 13 
or M 13 , since M 3 = n 2 M s /n 4 , and therefore M^M/M/ = n 2 2 M ls r M^~ r /n^ which makes 
every term in AiM/M/ a quartic in M 13> because M 3 is a linear function of M 13 and 
4— r is here zero or positive. 

(47) then shows that A 13 will be a quartic function of M 13 , if A x and A 3 are quartic 
functions of M x and M 3 respectively. But we know this to be the case for a single 
refracting surface. Hence it holds good of any system compounded of such surfaces. 

Now consider (48). If B 3 is a cubic in M 3 it becomes a cubic in M 13 . 

Again, if A 1 — M^ = a cubic ~U 1 in M^ 

M/(B 1 + A 1 M 3 / 1 // 8 ) = M/(B 1 {l+M 13 / 1 // 3 }+U 1 M 3 / 1 // 3 ) 

= M/M,/, (B^+U,//,), 

VOL. CCXXI. — -A. I 
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using (52), and by reasoning similar to the one given for A 13 the last expression is a 
cubic in M 13 . 

Hence B 13 is a cubic in M 13 . 

Consider now A 13 — M 13 B 13 . This is found, after some reductions, and using (52), 
to be 

(/a// 13 ) (A 3 -M 3 B 3 )+/3B 3 // 1 +M3M/(A 1 -M 1 B 1 ). 

Of the above terms, A 3 — M 3 B 3 is a cubic in M 3 and therefore also a cubic in M 13 . B 3 
is a cubic in M 13 . Aj— MjBj is a cubic in M l and when multiplied by M 3 M 3 2 becomes a 
cubic in M 13 . 

Hence if the condition (I) holds good for the components, it also holds good for the 
resultant system. But we have seen that it holds for a single refracting system ; 
thus it holds for any combination. Also B will be a cubic in M for any system. 

As regards C, examination of (51), remembering that for a single surface C 3 and G x 
are quadratics in M 3 , M 2 respectively, leads immediately to the conclusion that C is a 
quadratic in M for any system. 

We now come to the coefficient E. Here the single refracting surface gives no 
precedent for E x and E 3 . Let us examine the other terms in E 13 . These can be 
written in the form ./AM/M/ (B 3 - dAJdM^ + SA^if^+f.M^^f.A^G^ 
and, using Ai—BjMi = XJ l ; A 3 — B 3 M 3 = U 3 , where U 1? U 3 are then cubics in M l5 M 3 
respectively, this is found to reduce to 

-/AM 3 2 M 3 MM 3 c*B 3 /(M 3 +<^ (56) 



Now 



AxM/Mg 8 = quartic in M 13 . 
M 3 <2B 3 /eM 3 + dU 3 /cM 3 = cubic in M 3 = cubic in M 13 . 
M 3 2 M 3 (U x +/ 3 B 1 // % 13 ) = cubic in M 3 = cubic in M 13 . 

A 3 = quartic in M 13 . 
M 3 2 Ci = quadratic in M 13 . 

Hence the three terms in (56) are of form 

(quartic) (cubic) 4- (quartic) (cubic) + (quartic) (quadratic), 

and this leads to a rational integral polynomial of degree 7 in M 13 . 

Further consideration, however, shows that it is of degree 7 only in appearance, for 
the terms which can lead to expressions of degree 7 in M 13 are clearly 

-/ 1 M/M 3 3 A 1 d!B 3 / c «M 3 +3A 3 / 1 M/M 3 (U 1+ / 3 B 1 // 13 ), 
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or, dropping the factor fj 

M/M 3 {3A 3 {U 1 +f a B 1 lf w )-M 9 'A 1 dBjdNL a } 

= M 3 2 M 3 {9B^ 9 {V 1 +f a BJf a ) + S\J t (U 1 +f^ 1 lf n ) 
-Mg'MjBj dBJdM-s-M^ dBJdMJ. 

The term 

3U3M3 2 M 3 (U 1 +/ 3 B 1 // 1 3) 

is clearly of the form cubic x cubic and the terms leading to expressions of 7th degree 

reduce to 

M/M 3 {M 3 U 1 (3B 8 -M,dB a /dM,) + MA (aB*fdfi*-M a dB 3 fdM a )}, 

and since dM in /dM n ~f s /f m this can be written 

M 3 2 M 3 {M^ (3B 3 -M3dB 3 /^M3) + (M 3 B 1 /3// 1 3) (3B 3 -M 13 dBJdM„)}. 

But since B 3 is a cubic, in either M 3 or M 13 , 3B 3 ~M 3 (£B 3 /c£M 3 and 3B 3 — M. l3 dB 3 /dM ls 
are both quadratics in M 3 or M 13 . 

The above expression therefore reduces to M 3 x sum of two quantities each of form ; 
cubic in M 13 x quadratic in M 13 , that is, to a sextic in M 13 . 

We see, therefore, that those terms in E 13 which do not involve E : or E 3 are a 
polynomial of sixth degree in M 13 or M 13 . 

It follows that for a lens E is necessarily a sextic in the magnification. 

Suppose now that E x and E 3 are both sextics in M 1? M 3 respectively. Then E 3 is a 
sextic in M 13 and E 1 M 3 4 M 3 a will also be a sextic in M 13 , that is E 13 will again be a sextic 

• __ __ 
111 JXLjo* 

Hence, since any system is built up of combinations of lenses or single refracting 
surfaces, we find that E is a sextic polynomial in M for any system. 

Examination of particular cases shows that E is not, in general, divisible by A, so 
that the vanishing of the latter does not usually involve the disappearance of the 
second order terms. 

§9. Invariant Relations. 

Certain relations exist between the coefficients A, B, C, E which remain the same 

in form, whatever the number of refracting surfaces. One of these we have already 

dealt with, namely the fact that 

A-MB 

reduces to an expression of the third degree, i.e., the coefficients of highest degree in 
M in A and B are the same. 

This we shall refer to as the first invariant relation (I.). 

A second invariant relation takes the form 

B-C = £(l-MVHA/dM, . ....... (II) 

i 2 
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when we use tan /3 2n as argument, and 

B-C= -^{l-M 2 )+idA/dM, . (ir) 

when we use sin y 2n as argument. 

That these relations hold good for the single refracting surface is readily verified 

from equations (22), (23), (25), (27) and (32). Suppose now that, for systems 1 and 3 

separately, the relation 

B~C = o-(l-M 2 )+ic£A/(M 

holds, where <r = f or — ^ according to the nature of the argument, then from (48) 

and (51) 

B 13 -C 13 - B 3 -C3+M/(B 1 -C 1 ) + A 1 M 3 M3 2 / 1 // 3 

= (r (l-M 13 2 )^idA s /dM 3 +iM 3 2 (dA 1 /dM 1 + 4AM s fJf i )>- • (57) 
and from (47) 

f 13 dA 13 /dM 13 =/ 8 dAj<M ls +f 1 {M a 'K a a dASM 13 + A, d (M/M/)/dM 13 }. 

But 

dM 13 = (f 3 /f 13 ) dM 3 = (fjf 13 ) M 3 * dM x . 
Hence 

dA 13 /dM ls = dA.JdM 3 +M 3 * dAJdM, + (/, A//,) d (M 3 2 M 3 2 )/^M 3) 

and since M 3 /M 3 = const., the last differential coefficient is 4M 3 M 3 2 . 
Using this result (57) becomes 

B 13 -C 13 = <r(l-M 18 3 )+HA ia /dM 18 , 

which is of the same form as the equation we started from. Hence, if the two 
components of the compound system satisfy the second invariant relation, the 
resultant system also satisfies it. But we have seen that the relation holds good for 
single refracting surfaces — hence it holds good universally. 

It should be noted that the second invariant relation is really a first order relation 
and connects the first order aberration of the inclination of a ray, with the first order 
longitudinal spherical aberration. 

§10. The Constants A, B, C, Ejfor an Optical System Reversed and for Negative 

Lenses, 

Certain important general relations are found to hold between the constants 
A, B, C, E for rays going through an optical system and the corresponding constants 
A', B 7 , C, E 7 , for the same system reversed, and by making use of them we can obtain 
either set from the other. 

We arrive most simply at these relations as follows :— If after traversing the 
system we retrace our steps, the result is equivalent to compounding the system with 



SPHERICAL ABERRATION FOR A SYMMETRICAL OPTICAL SYSTEM. 57 

itself reversed, with the difference that, in the second set of refractions, the measure- 
ment of length parallel to the axis is reversed in direction. An examination of the 
equations (41) et seq.> § 7, on which the formulae of combination are based, shows that 
this is analytically equivalent to changes in the sign of the focal length in the second 
set of refractions. 
We have therefore 

/:=/, fs=-f, M 1 = M, M 3 =l/M, M 1 = M, M 3 =l/M, M ]3 = M 18 = 1, 

and we also find that f ls = oo. But /i 3 A 13 , and /i 3 E 13 have definite limiting values, 
and as jf 13 does not otherwise explicitly enter into the equations of combination, no 
difficulty arises on that account. 

Now, after retracing our steps in this way, we necessarily arrive at a perfect image, 
so that Ax± = and tan a 4 = tan /3 4 , leading to 

/i 3 E 13 = 0, 
and 

Bi 3 — C 13 = 0. 

These lead to the following identical relations 

A(M)/MT1P-A'(M- 1 ) = 0, (58) 

E {M)/WW-~E f (M- 1 ) + A (M) M~ 2 {B' (M" 1 ) M" 2 +3A / (M" 1 ) M^ + dA' (M" l )/dM} 

-A r (M- 1 )M- 2 {3B(M)-2C(M)} = 0, (59) 

B^M-VC^M-O+M-^ . . (60) 

Equation (58) may be written in either of the two forms 

A (M) = M 2 M 2 A' (M- 1 ) 
A(M)/<~M 4 A / (M- 1 )/n 2 2 . (Ill) 

This we shall refer to as the third invariant relation. It shows that, if we divide 

A by the square of the initial refractive index, the coefficients of powers of M 

equidistant from the beginning and end of the development are interchanged by 

reversing the system. Equation (59) becomes on multiplying up by M 2 M 4 , using (58) 

and simplifying 

E-M 2 M 4 E'+A{M 2 B , ~A/M-3B + 2C + dA/dM} EE0, . . . (61) 

omitting the arguments M, l/M of A, B, B', &c, since no confusion can occur. 
Use now the second invariant relation 

dk/dM = 4B-4C-4cr(l-M 2 ), 
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(6l) becomes 

E-M 2 M 4 E / + A{B + M 2 B / -2C~A/M-4(t(1~M 2 )} .= 0. , . . (62) 

Now substitute from (60) for B-f-M 2 B' the value 

C+IFC' + A/M, 

and (61) leads to 

E-M 2 M*E' + A{M 2 C'-C~4<r(l-M 2 )}EE0 ..... (63) 

For a single refracting surface, where E, E' are identically zero, this must lead to 

M 2 C'-C = 4<r(l-M 2 ) (IV) 

a result which is easily verified from equation (32). 

Now consider a system compounded of two systems. For the system direct, we 

have 

C 13 = C 3 + M 3 Oi. 

Similarly, for the system reversed, change C 3 into d, Cx into C ; 3 , M 3 into l/Mj. 

C 31 = C ! + Mx~ C 3 
M 13 2 C 31 - C 13 = M x 2 M/C, - M/d + M 3 2 C 3 - C 8 

= M/(M 1 2 C' 1 -C 1 )+M 3 2 C'3-C 3 

and using (IV) which we know to be true for a single refracting surface 

M 13 2 C 31 -C 13 = 4. [M 3 2 (1 _M X 2 ) + 1 _M 3 2 ] 

= 4o-(l~M 13 2 ). 

In other words (IV) will hold for the resultant system if it holds for the components, 
and therefore as in previous similar cases, it holds for any system. 

We shall call (IV) the fourth invariant relation. 

Equation (63) then shows that there exists a fifth invariant relation 

E = M 2 M 4 E ; . . (V) 

or 

E (M)/< = M 6 E ; (M~ l )/n 2 \ 

so that E possesses a property similar to that of A, previously noticed, viz., if we 
divide it by the square of the initial refractive index, the coefficients of powers of M 
equidistant from the beginning and end of the development in M are interchanged. 

Equation (59.) has therefore led us to two independent invariant relations. 

On the other hand it will be found that (60.) leads to no new relation. For if we 
substitute into it for B r — C and B — C in virtue of the second invariant relation, and 
then use the third relation, it becomes an identity. 
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It gives, however, on using (IV) to eliminate C, 



M 2 B' + B = A/M + 2C + 4<r ( 1 -M 2 ) 
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(64) 



which is a convenient form for calculating B'. 

If now A, B, C, E are known for any system, the corresponding quantities are 
immediately obtainable for the reversed system, A', B', C, E' being given by 
equations (III), (64), (IV) and (V) respectively. 

This will generally halve the labour of calculations, if it is found desirable to 
tabulate these constants for a complete set of lenses. It will then be sufficient to 
start from the equi- convex lens and vary the curvatures in one sense only. 

Incidentally we note also that the second invariant relation enables us to find C, 
so soon as A and B are known, so that only A, B and E require to be calculated. 

The aberration constants for a reversed system have a further important application 
in the case of lenses. Consider a positive lens (fig. 4), the initial ray converging to 
I and the final ray to I 4 . If now we interchange the full and dotted portions of the 
initial and final rays in fig. 4, we obtain, since here the initial and final media are the 




Fig. 4. 



same, the case of a ray going through a lens in which the front and back character of 
the two surfaces have been interchanged. In fact n and r 3 have been interchanged 
and the sign of the thickness c 2 has been reversed. This leads to a negative lens, of 
the same numerical power as the original positive lens, and with the same mean 
curvature, but a negative thickness. Such a lens, of course, is not physically 
realisable, although a part of it can be physically obtained by rotating the wedge 
beyond the intersection V of the two surfaces. 
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But, in the case of the ideally thin lenses, where the thickness is zero, the ideally 
thin positive and negative lenses, having the same mean curvature and numerical 
power, correspond in this way. 

Now I 4 and I are also interchanged. If we consider I 4 as the initial point, then we 
are really considering a set of rays starting from I 4 in the last medium and travelling 
backwards through the original positive lens. In other words, the aberration constants 
for the corresponding negative lens are identical with those for the original positive 
lens reversed, and the equations (III), (64), (I\^) and (V) are applicable to calculate 
them. 

This, again, will greatly diminish the work of calculation. In the case of ideally 
thin negative lenses, we see that A, B, C, E are directly obtained from the corre- 
sponding thin positive lenses. In the case of thick negative lenses the corresponding 
positive lens has a negative thickness. 

Now for various reasons it will probably be convenient, in calculating A, B, C, E 
for lenses, to express them in the form 



A c = A + c (dA/dc) 



0> 



&c, where c, the thickness, is small, as it usually is in practice, and A refers to an 
ideally thin lens. 

When the formulae are put in this form, it is perfectly simple to calculate A_ c , B_ c , 
&c, and then to obtain the corresponding results for the negative lens with a positive 
thickness. 

§ 11. Explicit Values of A., B, 0, E/or a Thick Lens (Tangent Formula). 

For the purposes of numerical calculation and comparison with correct trigonometri- 
cally found values, we have worked out explicitly the form of the expressions A, B, 
C, E for a thick lens, when we use tan/3 4 as the argument; the formula are 
expressed in terms of the focal lengths of each surface and of the combination and the 
thickness does not appear explicitly. The initial and final media being the same 
n =n 4 and we have written n = n 2 /n . 

The work of algebraic calculation has been straightforward but extremely heavy, 
and we therefore omit it here entirely, the object being to publish the results for 
reference, in case other workers desire to use them for tabulation purposes, but it is 
hardly to be expected that designing opticians should work direct from the algebraic 
expressions as they stand. 

For the purpose of this section we shall write 

A = A + AiM + A 2 M 2 + A 3 M 8 + A 4 M 4 . 
B = B + B 3 M + B 2 M 2 + B 3 M 3 . 

C-C + C 1 M + C 2 M 2 . 

E = E +E 1 M + E 2 M 2 +E 3 M 3 +E 4 M 4 +E 5 M 5 + E 6 M tt . 
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The suffixes here have a different meaning to that which has been previously 
ascribed to them, but no confusion is likely to arise on this account. 

The values are as follows, f denoting the focal length of the lens, f ly f 3 the focal 
lengths of the surfaces, so that 



An — 



-0 



A 1 = 



2,71? (n— 1) 
1 



/i = n/in-1), f 3 = - r J(n-l) 



A 2 = 



B„ = 



2nMn-i ) 2 r 1 ^^ /3+4(w+1)2 ff~ 4w(1+w+wii) i +w2(w+1)a 

2re 2 (w-l) 2 L / 3 2 / 3 2 7 /i/ 3 \/i /s'J 

Ai with/i and / 3 interchanged 

A 

1 



'0 






B 1 = 

B 2 = 

B 3 = 

O = 

c 2 = 



6/Vl , 1\ 



/ 



,./l 1 
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"C 1 ___ 
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f 

n 



2n ift + ft) + A{n+lY (l 2+ k°) + 

2n(n+iy(£+£j + [_(n+iy 



12 (n+l) 2 {n a + n+ 1} 
J1J3 



+ 4:(n 2 + n + l)(n* + 3n+l)](j+j)\ 



+ 4,n(n 2 + n + l) (4+^] + 2 {n + l) 2 (n 2 + 3n+l) 

\/ 3 fJ 



2n 2 ( n ^if(M£) + 2n"^ 

f J 



E 4 = E 2 with f x and / 3 interchanged. 



5 — -^1 J? 



?) 



>> 



>? 



6 — -^0 J 



>J ?> 



J ? 



?> 
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It should be noted carefully that all the above refer to expressions in terms of the 
tangent of the Gaussian inclination, this being the argument we have used in the 
numerical work. 

§ 12. Values of A, B, C, TLfor a Thin Lens. 

When the lens is thin, we have the relation 

J: -+.! = ! 
J i Ja J 

which enables the values of § 1 1 to be considerably simplified. 
In this case it is useful to introduce a quantity K such that 






mean curvature of the lens _ f /1_ 1 
power of the lens 2 \r 1 r 3/ 



When this is done the constants A, B, C, E take the following forms : — 

A= ~{{l-M) 2 /2n} {(» + 2)[(l-M)K-(l+M)(w + l)/(« + 2)] 2 

+ « 3 (l-M) 2 /4(w-l) 2 -w 2 (l+M) 2 /4(n+2)}. 

B = (l-M) 2 K 2 {«-l-(w + 2)M}/2w + (l-M)K(l + M + 4M 2 )(n+l)/4» 

+ (l-M){M(l+M)/4» + (l-M)[3»-2n 2 -3 + M(6w--4w 2 -3)]/8(n--l) 2 }. 

C= -3(l-M) 2 K72w + 3(l-M 2 )K(n+l)/4«-f(l-M 2 )-fn(l-M)7(«-l) 2 . 



E 



3(1-M) 2 

128(w-l) 4 » 3 



"(l + M.' 4 ) (-4ri 6 + 8w 4 -n 3 -4n 2 +3n-l) 

+ M(l + M 2 )(8w 6 -16n 5 + 4n 4 + 4w 3 -12w 2 +12n-4) 
+ M 2 (-lGra 7 + 32n 6 -8w 5 -8n 4 +10n 3 -16n 2 +18w-6) 
+ 8 (n-l) (1-M) (1+M) K {(1 + M) 2 

(-2n b + 5n 4 ~2n 3 -3n 2 + 3n-l) 
+ M (2w 6 -8w* + 4» 3 + 2w 2 )} 
+ 8(«-l) 2 (l~M) 2 K 2 {(l + M) 2 

(-2n 6 + 8» 4 -7« 3 -6w 2 + 9«-3) 
-2Mw 2 (w 2 -2n-l)} 
+ 16(»-l) 3 (l+M)(l-M) 3 K 3 {2w 4 -4w s -2n 2 + 6»-2} 
+ 16 (w-l) 4 (1-M) 4 K 4 (-n 3 +3w-l) 

We notice that when M = 1 (which gives one of the zeros of A) B, C and E all 
vanish with it, and also E/A remains finite. Hence in this case, the term E£ 4 4 will 
not rise in importance, even when A = 0. But in this case A may have two other 
real zeros, and these are not zeros of E, so that E plays an important part in the 
neighbourhood of such zeros. 

k 2 
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§ 13. Numerical Test for a Single Lens. 

To test the formulae, a number of longitudinal aberrations were calculated trigono- 
metrically for five positive lenses of refractive index 1*52, unit focal length and 
thickness ^. The first was a meniscus-shaped lens for which r 3 = 0*349418, r 3 = 1. 

The second was a plano-convex lens, of which the convex side is towards the 
incoming light. The third was an equi-convex lens. The fourth and fifth were the 
second and first reversed. If K has the meaning defined in § 12, the values of K for 
these five lenses are 1*93094, 0*96154, 0, -0*96154 and - 1*93094 respectively, so 
that they proceed by approximately equal steps of K. 

The constants A, B, E were calculated for these five lenses and the longitudinal 
aberrations computed from the formula. The rays selected met the first principal 
plane at a distance 0*15 from the axis, corresponding to an apertures/3 '4, nearly. 

The results are shown in Table III. 



Table III. — Longitudinal Aberrations of Five Selected Lenses. 



M. 


Lens 1. 


Lens 2. 


Lens 3. 


Formula. 


True. 


Per- 
centage 
error. 


Formula. 


True. 


Per- 
centage 
error. 


Formula. 


True. 


Per- 
centage 
error. 


3 

2 

0-5 


-0-5 
-1 
-2 


-7- 

-1- 
-0- 

-0 
-0 
-0 

-1 


98683 
55287 
003527 
•06028 
•20902 
•42692 
•00383 


-7-87422 

-1-52334 

-0-003530 

- • 06043 

-0-20946 

-0-42772 

-1-00531 


1-41 
1-94 
0-09 
0-24 
0-21 
0-19 
0-15 


-1-12450 

-0-40167 

-0-011697 

-0-024728 

-0-083863 

-0-18428 

-0-48874 


-1-12814 

-0-40283 

-0-011705 

-0-024749 

-0-083950 

-0-18448 

-0-48930 


0-32 
0-29 
0-07 
0-08 

o-io 

0-11 
0-12 


-0-32378 

-0-14240 

-0-033008 

-0-036194 

-0*057067 

-0-094736 

-0-21707 


-0-32439 

-0*14282 

-0-033096 

-0-036213 

-0-057101 

-0-094787 

-0*21724 


0-19 
0-29 
0-27 
0-05 
0-06 
0-05 
0-08 



M. 


Lens 4. 


Lens 5. 


Formula, 


True. . . 


Percentage 
error. 


Formula. 


True. 


Percentage 
error. 


3 

2 

0-5 


-0-5 
. -1 
-2 


-0-09S969 

-0-046089 

-0-073959 

-0-099843 

-0-13967 

-0-19084 

-0-32454 


-0' 

-0' 

-0 

-0 

-0 

-0 

-0 


099086 
046128 
074535 
•100044 
•13984 
•19101 
•32476 


o- 

0' 
0' 


. 




12 
30 

77 
•20 

12 
•09 
•07 


-0-10401 

-0-013671 

-0-14762 

-0-22172 

-0-32828 

-0-45753 

-0-77417 


-0-10444 

-0-013686 

-0-15068 

-0-22324 

-0-32966 

-0-45887 

-0-77511 


o- 

0' 

2 

0' 






41 

11 

03 

68 

•42 

•29 

•12 



The mean percentage error of these results is 0*34, so that the formula determines, 
on the average, the longitudinal aberration correct to 1 part in 300. The bulk of 
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this, however, is contributed by three cases, namely : M = 3 and 2 for lens 1 where 
the aberrations are very large and differ very widely from the usual first and second 
order approximations, so that, although the error of the formula approaches 2 per cent., 
it nevertheless represents a great improvement upon these approximations ; and 
M = 0*5 for lens 5, which corresponds to extreme curvature and highest inclination, 
so that one of the angles of refraction is as great as 48^* degrees. Even here the 
table below shows that the formula is an appreciable improvement on the usual 
second-order approximation. If these three cases are omitted, the mean percentage 
error works out to be about 0*21, so that in general the formula determines the 
longitudinal aberration correct to about 1 part in 500. 

It is interesting to note what the usual first and second order approximations lead 
to in a few cases. 









First order 


Percentage 


Second order 


Percentage 


Lens 


1, M = 


3 . . . 


approximation. 


error. 


approximation. 


error. 


-1-34012 


83-0 


-2-45432 


68-8 


5J 


1, M = 


A ... 


-0-48545 


68-1 


-0-82104 


46*1 


)> 


1, M = 


~~ 2i . . . 


-1-348142 


34-1 


-0-88456 


12*0 


>) 


2,M = 


o ... 


-0-67356 


40-3 


-0-94436 


16-3 


»» 


5, M = 


0*5. . . 


-0-12136 


19-5 


-0-14307 


5*1 



This gives a measure of the numerical improvement effected by the fractional formula 
whenever the usual method of approximation is seriously out, even though in none of 
the cases above does the convergency of the series actually fail. 

In the above the series are in powers of tan (3 2 . Had they been taken in powers of 
tan a 2 , as is frequently done, the first and second order approximations would have 
been far worse. 

One interesting outcome of these calculations relates to the relative importance of 
the terms in E£ 4 4 and At*. The ratio E£//A is small in every case taken (of course 
these exclude the neighbourhood of points where A = 0, where naturally E becomes 
of great importance). But for the set of magnifications taken, the greatest ratio of 
the second term to the first is less than 0*03 and the mean value of this ratio is only 
0*0082, so that, in fact, the E term — although so complicated algebraically — does not 
exercise any great influence numerically. 

This is important, as it shows that, at any rate for lenses, it does not require to be 
computed with anything like the same order of accuracy which is needed for A and B. 



§ 14. The Singular Inclination and Convergency Factor for any System. 

Referring again to fig. 2 we see that A = a 2 and F I = — Ax 2 for rays proceeding 
through the system reversed and initially parallel. 

Thus using accents, as before, to denote the coefficients and inclinations for 
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the system reversed, and noting that the accented coefficients all refer to zero 
magnification, we have, using tangents 



and 



tan X = tan /3' 2 ( 1 + B' tan 2 p 2 )f( 1 + C tan 2 /3' 2 ) 
-F I =fn (A' tan 2 /3' 2 + E' tan 4 /3' 2 )/(l + B' tan 2 /3' 2 ). 



Here the suffixes in the A, B, C, &c., A 7 , B 7 , C, &c., have the same meaning as 

in §11. 

Thus 

-1/M = (A' tan 2 /3' 2 + E' tan 4 /3' 2 )/(l + B' tan 2 /3' 2 ), 

whence, developing cot 2 /3' 2 in descending powers of M and stopping at the second 

term 

cot 2 /3' 2 = -A' M-B' + E' /A' . 

Substituting into 

cot 2 X = cot 2 /3' 2 + 2 (C'o-B'o) 

which is valid to the same order of approximation, we obtain 

cot 2 X= -A / M + 2(C / -B / )-~B , +EVA / (67) 

as the second approximation for the singular inclination when M is large, the first 
approximation being cot 2 X = — A' M. 

To the same order the convergency factor is 

1 + tan 2 a (A / M+ {2 (B' -C ) +B' - (E' /A' )}), 



('. o. . 



l4-(n 2 2 tan 2 /3 2 K 2 ){A / M 3 +(3B' -2C' -E / n /A , )M 2 }, . . . . (68) 



/3 2 referring to a ray passing through the system in the standard sense. 

Now, using the equations (III), (64), (IV) and (V) of § 10 and equating suitable 
coefficients, we find that 



E' = K7»/) E 6) 

B'o = (», 2 K) (A 3 +2C 2 -B 2 ) -4<r 



A 4 — B 3? 



A o — (^0 l n 2 ) A 4 , 

C' = (n 2 /n 2 *)a-^ 

3A 3 = 4B 2 -4C 2 + 4«K 2 )^ 

whence, after substitution, (68) becomes 

1 + tan 2 fa (B 3 M 3 + {B 2 -E 6 /A 4 } M 2 ). ....... (69) 

Now, if our B leads to a sound approximation to the convergency factor for M large, 
this should be 

1 + B tan 2 j3 2 , 
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or, to the same approximation which we have been using 

l+tan 2 /MB 3 M 3 + B 2 M 2 ) (70) 

We see, therefore, that the development of the correct convergency factor in descending 
powers of M will give a result which always agrees with our B, so far as the highest 
term in M is concerned, but makes the term in M 2 in general different. 

In the case of a lens E 6 /A 4 is in general small, compared with B 2 , so that this 
discrepancy makes little difference, but it may well be that, when we come to deal 
with more complicated systems, this will not be the case. 

A little consideration, however, shows that, when this is so, our formula is very 
readily corrected so as to take this difficulty into account, without involving any 
lengthy numerical computation. 

If we consider the formula 

Ax 2 = n 2 /{A* 2 2 + (E-AE 6 M 2 /A 4 )* 2 4 }/{1 + (B-E 6 M 2 /A 4 )* 2 2 } 

it is clear that it leaves the development of Ax 2 in powers of t 2 unaltered as far as 
the second order inclusive. It alters the coefficient B 2 of B so as to make the two 
leading terms agree with (69). It also alters E in such a way as to remove the term 
in M 6 and reduce E to a quintic. In fact it gives for the new E the remainder obtained 
after the first step in the division of E by A, according to the usual process. 
In practice, the terms in (E 6 /A 4 ) M 2 are very readily added as follows : — 

Ax _ „ , M 2 + E^ /( 1^E 6 M% 2 /A 4 ) 
2 - n 2 j f+fyf/fr^wtf/j^ 

and this amounts to applying the same corrective factor l/(l— E 6 M 2 £//A 4 ) or 
l/(l—'En % 2 /n 2 2 A 4 ) to the second terms in both numerator and denominator. This 
factor, expressed in terms of the inclination of the incident ray, is independent of the 
magnification, and a short table will enable it to be found in any given case without 
difficulty. 

A similar correction has then to be made in C ; in order to keep the development 
of tan a 2 the same we must have 

tan a 2 = t, { 1 + (B-E G M 2 /A 4 ) t>} / { 1 + (C-E 6 M 2 /A 4 ) t>}, 

and writing this as 

^{1 + B^ 2 2 /(1~E 6 M% 2 /A 4 )} 
l + a 2 2 /(l~E 6 MV/A 4 ) 

we see that the same corrective factor has to be applied to all the second terms in the 
formulae. 

In the above we have used tan /3 2 as our argument, but the formulae and the 
correction take precisely the same form if sin y 2 is the argument. 
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If we apply this corrective factor to the first two entries of Table III., which give 
a large percentage error — -these correspond to cases approaching the failure of 
convergency and are therefore critical, we find, for lens ( I ) 



M. 


True 
aberration. 

-7-87422 
-1-52334 


Formula 
uncorrected. 


Percentage 
error. 


Formula 
corrected. 


Percentage 
error. 


3 

2 


-7-98683 
-■ 1 ' 55287 


1-41 
1-94 


-7-90811 
-1-54089 


0-43 
1-16 



which shows a very sensible improvement. 

The significance of this alteration is brought out more clearly when we consider 
the limiting case M = o° that is, rays actually issuing from the front focus (the 
case of an eye-piece). In this case, the geometrical image being at infinity, it is 
inconvenient to define the emergent ray by means of either longitudinal or transverse 
aberrations. 

Let us consider the intercept of the ray on the back focal plane. 



Also 



where 



This = (—n 2 /M + Ax 2 ) tan a 2 

= n 2 f[ (A tan 2 /3 2 + e tan 4 /3 2 )/(l + b tan 2 /3 2 ) — M] tan a 2 . 

tan a 2 = (tan a /M) (l + bM~ 2 tan 2 a )/(l +cK~ 2 tan 2 a ), 

b = B-E 6 M 2 /A 4 , c = C-EeM7A 4> e = E~E 6 AMVA 4 , 



and it is clear that, in the limit, where M (and M) = oo, tan <x 2 must be finite. 

This requires that b and c shall be of order M 3 and M 2 respectively, which is right, 

and leads to 

tan a 2 = nb s tan 3 <x fn ( 1 + c 2 tan 2 a ), 

b s and c 2 being the coefficients of M 3 and M 2 in b and c respectively. 

On the other hand, it is equally obvious that the intercept on the back focal plane 
must also approach a definite limit. Hence the factor 



i.e., 



or 



(A tan 2 p 2 + e tan 4 &)/( 1 + b tan 2 fa) - M, 

- 1 + (A-M5)tan 2 &/M + 6 tan 4 /3 2 /M 
l/M + 6tan 2 /3 2 /M 

- 1 + ( A - Mb) tan 2 a n /MM 2 + e tan 4 a /MM 4 

l/M + 6tan a a /MM a 

must tend to a finite limit as M approaches co . 

This necessarily involves that (l) A-M6 has M 3 in its leading term — a result 
already established, and (2) that e involves M 5 (and not M 6 ) in its leading term. 
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Hence, whenever we deal with incident rays actually passing through the front 
focus — and this necessarily occurs as soon as the results of the present paper are 
applied to aberrations off the axis — the modified B, C and E have to be used. 

§ 15. Combination Formulas for More than Two Systems. 

The combination formulae (47), (48), (51) are capable of explicit generalisation for 
any number of systems. 

In the case of A and C successive applications of equations (47) and (51) lead at 
once to the results 

«/135...2n + l"»135...2n + l == (/l^-l) M. 35 ... 2w + 1 JM. 35 ,. <2n + 1 

+ VJ3A3) -M. 5...2n + l-M- 5...2n + l+ ••• 

+ \t/2n -lA-2n-l) -M- 2n+l«™- 2n + l "h/2n + l^-2n + l« • \'l) 

^135...2n + l ==: ^1-™- 35..-2n + l + W"- 5...2n+l"+" ••• 

+ C 2M _iM 3 2n+1 + C 2n+ i. . . . . . . . (72) 

in which we have reverted to the notation of §§ 2, 7. 
In the case of B, we have for three systems 



B 135 = B^+B&^+AM^A/fu 

= B B +B 8 M 5 » + B 1 M 86 8 +(/ 1 A 1 ) Ma^^+C/A) M 5 M/// £ 



35 

r TWT 2//* _!_/./• A MV/TM-2// 

5? 



and the general law of combination is at once obvious. 
We have 

Bi3 5 ...2n+i — BiMgg 2n+1 + B 3 M 5 2n+1 + ... + B^^M 2n+1 + B 2w+1 

+(f 1 A 1 ) (MMV/) 35 ... 2 „ +1 +(/,A 3 )(MMV/),., n+1+ ... +(A_ 1 A 2 ,_ 1 ) (MM 2 //k +1 - (73) 

The case of E is more difficult and we have not been able to obtain a form in which 
it can be written down for the combination of n systems. 

But we can deal with it as follows, by determining the contribution of any one 
component to the whole : — 

Consider three systems 1, 3, 5. It will be convenient to look upon 3 as a single 
lens, forming part of a larger system. 1 will then be the system of lenses preceding 3, 
and 5 the system following. 

Applying equation (49), simplified by the use of the second invariant relation, we 
find 

/ 135 E 135 =/ 6 E 6 +/ I3 E 13 M 5 2 M^ 

+/ 13 A 13 M/M 6 2 (3A 6 /M 5 -3B 5 +4C 5 +4o-{l-M/}) 

+ 3B 13 / 5 A 5 M 5 3 -2C 13 / 5 A 5 M 6 2 . 

VOL. CCXXI.— A. L 
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Applying the equations of combination a second time, and picking out the terms 
involving A 3 , B 3 , C 3 , E 3 , we find these to be 

/3E3M/M/+M/M/ [/AM/M/ (3A 3 /M 3 -3B 3 +4C 3 )+/ 3 A 3 M 3 2 (36,-20,)] 
+/ 3 A 3 M 5 2 M 5 2 [3A 6 /M 5 -3B 5 + 4C 5 + 4<r (l -M 6 2 )] 
+ 3B 3 / 6 A 6 M ft 2 -2C 3 / 6 A 5 M/. 

Hence the lens 3 contributes to the final E 



2 
L 5 



/, A s { 3M 5 M/ [ A 6 + (fjf 3 ) A.M^M,/] + 4<r ( 1 - M 5 2 ) M/M E 

+ M/M/[-3B 5 + 4C 5 + (3B 1 -2C 1 )My} 

+ B 8 M 6 8 (3/A-3/iA 1 M 36 2 M 8B 8 ) 

+ C 3 M/(-2AA,+ 4/ 1 A 1 M 35 2 M 38 2 ) 
+/ 8 E 8 M 5 *M ( 



4 

u 5 ? 



and the A's, B's and C's in the curled brackets can be expressed in terms of the 
individual lenses of the system by means of equations (71), (72), (73). 

If we denote the coefficients of A 3 , B 3 , C 3 , E 3 , in the above by Z 3 , m 3 , _p 3 , g 3 , then the 
contribution of the individual lens to E 

= Z 3 A 3 +m 3 B 3 +jp 3 C 3 + g 3 E 3 . 

Hence, if we vary K 3 for this lens, keeping focal length and magnifications 
unaltered 

AE = AK 3 (Z 3 3A 3 /3K 3 +m 3 3B 3 /3K 3 +p 3 3C 3 /3K 3 + g 3 3E 3 /3K 3 ). 

If all the lenses are simultaneously varied, then we have 

AE = 2 AK (1 3 A/3K + m 3B/3K +p 3C/3K + q 3E/3K). 

"We have similar equations for AA, AB, AC, but they take a much simpler form. 

Using these, we can, if we have enough lenses, vary the K's so that, between limits, 
we can make our four constants A, B, C, E take up any assigned values, or, if we 
wish to keep any one constant whilst slightly varying the others, we have a linear 
relation between the AK's. 

§ 16. Conclusion. 

We have now established a formula of fractional type for the longitudinal aberration 
of a symmetrical system which, while algebraically correct as far as the second order, 
does in fact, give results beyond this order in those numerical cases which have been 
tried, and largely overcomes the difficulties of slow convergency in critical regions. 
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We have further obtained a method for calculating the coefficients of this formula 
for any symmetrical optical system in terms of the coefficients for the components > 
in such a way that the effect of any single component upon the whole combination is 
immediately obtained. 

In considering the convergency of the series usually employed, we have found that 
the value of the approximation depends upon the particular variable employed, and 
that if we wish to avoid trouble owing to lack of convergency we must use sin a or 
tan a (or a suitable multiple of these) as argument, where a is the inclination of the 
original incident ray. 

The numerical success of the new formula appears to suggest that progress in the 
algebraic treatment of symmetrical instruments is to be sought, not so much along 
the lines of developments in series, but in other mathematical directions such as 
continued products, or possibly continued fractions. 

The next step would be to develop the 'method so as to cover the second order 
approximations to the emergent inclination. This will enable us to deal with aberra- 
tions off the axis of the system. 

Some progress has already been made by the authors in this direction . and the 
results will, it is hoped, form the subject of a later communication. 



